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Abstract 



We describe some new families of quasimodes for the Laplacian perturbed by the 
addition of a potential formally described by a Dirac delta function. As an application 
we find, under some additional hypotheses on the spectrum, subsequences of eigen- 
functions of Seba billiards that localise around a pair of unperturbed eigenfunctions. 

1 Introduction 

One of the unsolved questions in the analysis of quantum eigenfunctions concerns pos- 
sible limiting distributions as the eigenvalue tends to infinity. For eigenfunctions of the 
Laplace operator on certain surfaces with arithmetical properties, it has been proved [U [2] 
that there is only one possible limit; all sequences of eigenfunctions become uniformly 
distributed. On the other hand, Hassell [31 has proved the existence of chaotic billiard 
domains in M? for which zero-density subsequences of eigenfunctions fail to equi-distribute 
in the limit. 

We consider the Laplace operator plus potential supported at a single point. Such 
a potential has been variously referred to as delta-interaction potential, Fermi pseudo- 
potential or zero-range potential in different parts of the literature. Mathematically this 
operator can be constructed using the tools of self-adjoint extension theory. 

We will prove our results for the case where the underlying space is a compact 2- 
dimensional manifold, for which the Laplace operator has eigenfunctions and eigenvalues 
denoted by <Pj and Ej respectively. We perturb this operator with a delta potential sup- 
ported at the point p, which will remain fixed throughout, and suppressed from notations. 
This perturbation can be realised by a 1-parameter family of self-adjoint operators, H@, 
indexed by an angle which controls the strength of the perturbation. 

We fix a finite interval ICR containing at least one Ej, and define, for notational 
convenience, 





Let a 6 [0,1]. We define 



for /i a solution to 

oa.,)^E(^-r^le)?ff- < 3 » 

Our main results are as follows: 
Theorem 1.1. The pair (if), ft) is a quasimode for H@ with discrepancy d, where 

(1 - a) 2 C/(0, /*) + ^ 2 E %0/ (l + ^A* + " M)V SSft 
d 2 = ^ g 7 1 jJ (4) 

Cl(2,M) + ^E^(^-i^) 



Furthermore, iftpx,tp2 are defined by ^ for /ii 7^ two solutions of Q th 



en 



The construction of families of quasimodes is a key step in Hassell's proof [3], as 
well as the proofs of many recent results on localisation of quantum eigenfunctions [U 
13 |9] . One reason for this is that quasimodes can often be used to approximate 
eigenfunctions. In general (see the introduction to section [3] for precise statements) the 
smaller the discrepancy, the closer quasimodes are to true eigenfunctions. For this reason 
it is important to know when the discrepancy can be made small. In this direction we 
have the following corollary to theorem 

Corollary 1.2. Let a = 1 and let I = [0, T] where T > E\. Then the discrepancy d of 
the quasimode ip satisfies 

d«^=. (6) 



Let (7 = and I be any interval containing at least two Ej. If fx £ I then the discrepancy 
d of ip satisfies 

d^^t(I), (7) 
where £(I) is the length of I. If, additionally, I contains precisely two Ej, then we have 

d<\m (8) 

In particular, the quasimodes with a = 1 and [i held fixed or slowly growing, can be 
made arbitrarily precise by choosing T as large as desired. 

We are interested in ascertaining when true eigenfunctions of Hq have mass supported 
on our quasimodes. Without any assumptions on the spectrum of the Laplacian we can 
prove the following. 
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Proposition 1.3. For any consecutive eigenvalues E a < E\> < E c < from the sequence 
{Ej)JL l; let I = [Eb,E c ] and take a = 0. Choose so that n G I. Then there is an 
eigenf unction (j) of Hq with eigenvalue in the interval (Eg,, E&) such that 

1/2 



Proposition 1.3 is most interesting when the sequence of eigenfunctions <Pj do not equi- 
distribute. (For example if they are solutions to a PDE which is subject to separation of 
variables; see below.) Then, by considering an infinite subset of the spectrum {Ej} along 



which the right-hand side of (|9j) is bounded away from zero, proposition 1.3 proves the 
existence of a sequence of eigenfunctions of Hq which fail to equi-distribute. Such a subset 
of {Ej} does exist since the mean level spacing is constant. 

Clearly the best that proposition |1.3| can achieve is to prove that a sequence of quasi- 
modes has an overlap of up to l/\/3 with a subsequence of true eigenfunctions. In order 
to prove that a sequence of quasimodes converges fully towards a sequence of eigenvalues 
of Hq we need to make some assumptions on the spectrum of the Laplacian. Sufficient 



conditions for this and a precise statement of the result (theorem 4.4) are given in section 

H 

The history of the study of the spectral properties of differential operators perturbed 
by the addition of a delta scatterer goes back at least to [10J , in which a one-dimensional 
lattice of delta interactions was used to model an electron moving in a crystal lattice. A 
comprehensive historical review is given in the appendix to the book |llj . 

Part of our interest in the subject comes from the Seba billiard which was introduced 
in [12]. In this work, a hard- walled rectangular billiard with a potential supported at a 
single point was considered. In terms of classical dynamics, the motion is integrable, since 
only a zero-measure set of trajectories meet the point at which the potential is supported. 
However, diffraction effects are introduced when one considers the quantum spectrum of 
the corresponding Schrodinger operator. 

Seba billiards have become important since the observation |13l IT4"] that the quan- 
tum spectral statistics belong to a new universality class, different from the classes from 
random matrix theory conjecturally associated to chaotic dynamical systems |15} [T6] or 
the statistics of a Poisson process conjecturally associated to fully-integrable dynamical 
systems [17J. It is now known that general integrable systems perturbed by the addition 
of such a localised scatterer also belong to the same universality class [18], as do quan- 
tum Neumann star graphs [HH ED]. Characteristic features of the spectral statistics of 
this universality class are an exponential decay of large level spacings, together with level 
repulsion. 

Several analytical studies of these spectral statistics have been made [T51 |2T | 122 ] [23 } |24" 1 
125] . Typically, a key feature of these arguments is the assumption of Poissonian behaviour 
for the eigenvalues of the billiard table without scatterer, a conjectured consequence of 
the integrable dynamics (the Berry- Tabor conjecture) |17| . 



In the final section of this article we apply theorem 4.4 to the original Seba billiard. Our 
final result is a proof that there exists a subsequence of eigenfunctions of the Seba billiard 
that become localised on a pair of consecutive eigenfunctions of the unperturbed billiard, 
if the spectrum of the unperturbed billiard satisfies an assumption which is consistent with 
the Berry- Tabor conjecture. 
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This result is a rigorous derivation of a formal argument first proposed in [26] and 
mirrors a related result proved for quantum graphs with a star-shaped connectivity [27J. 
These so-called quantum star graph can be considered as a singular perturbation of a 
disconnected set of one-dimensional bonds, each supporting a wave-function. In |27j the 
existence of subsequences of eigenfunctions that become localised on a pair of bonds was 
proved. This is exactly analagous to the localisation onto a pair of unperturbed billiard 
eigenfunctions in theorem 4.4 In both |27| and theorem 4.4 the main idea of the proof is 
to show localisation in an eigenfunction with eigenvalue lying between two closely-spaced 
eigenvalues of the unperturbed problem. 



2 Realisation of the perturbed operator 

Let Ai be a compact 2-dimensional Riemannian manifold, possibly with piecewise-smooth 
boundary, and let A be a self-adjoint Laplacian on Ai. 

The realisation of the operator formally defined by 

H = -A + c5(x-p), (10) 

where p £ Ai and 5 is the Dirac delta function, using the theory of self-adjoint extensions 
is given in many places in the literature. We refer the reader to |28|, [29] for the details. 
Here we recapitulate only that which is necessary to fix notations. We denote by || • || and 
(•, •) the norm and inner product of L 2 (Ai). 

Since Ai is compact, —A has a complete basis of eigenfunctions, with corresponding 
eigenvalues Ej which we write in non-decreasing order. 

We will remove from the list of eigenvalues any Ej for which <Pj(p) = 0. Such eigen- 
functions are not affected by a delta-scatterer at p, and so it is convenient to exclude them 
from the spectrum. This further allows us to assume that the spectrum {Ej} is simple, 
without losing generality. 

To see this, consider an eigenspace of dimension r > 1 spanned by the eigenfunctions 
{01, . . -4>r}- Then the vectors ((f)i(p), ■ ■ ■ ,4> r (p)) T and (R, 0, . . . ,0) T in C r , where 

have identical norm. This means that we can find a unitary r x r matrix mapping the first 
vector to the second. Multiplying U by the vector of eigenfunctions (<p\, . . . ,<?V) T leads 
to a new basis for the eigenspace, in which all but the first eigenfunction vanishes at the 
point p, and the corresponding eigenvalue is counted with multiplicity one. 

The resulting spectrum is therefore ordered so that 

Ei<E 2 <E 3 --- (12) 

We will frequently use Weyl's law with remainder estimate |30| : 

N(E) := £ |^-(p)| 2 = ^ + O^ 1 / 2 ), (13) 
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where the implied constant]^] may depend on the position of the point p G M . 
Define 



oo 



9±i(l):=E ?ife)|M. (14) 

j=l 3 + 1 

Then g±\ S L 2 (A4) and in fact they are the Green functions for the resolvent of —A at 
the imaginary energies ±i, satisfying 

(/,< 7±i ) = (-A±i)- 1 /(p). (15) 

In particular, 



(# i ,g±i) = ((-A±i)- 1 #,)(p) 



£,±i' 



(16) 



which will be useful to know later. 
Let 

D p :={/eDom(A):(/,5 p ) = 0}, (17) 
and define the operator Hq with domain D p by 

H : / ^ -A/. (18) 

i^o is a symmetric, but not self-adjoint operator. In fact its deficiency subspaces are 
spanned by g±\. 

It follows from the von Neumann theory [31 j that, 



Dom(i/ *) = V p span{ 5i , g^}. (19) 

Since the deficiency indices are equal, Hq possesses self-adjoint extensions, constructed as 
follows. 

First of all, note that we can write for ifi G Doul^q), 

V> = ^ + a + (V)5i + MV03-i, (20) 
where if) G T> p and a±{ip) G C. In fact we have 

# *V> = H i> + ia + (^) gi - ia_(V>)S-i- (21) 

Since the deficiency indices of Hq are both equal to 1, there is a 1-parameter family of 
self-adjoint extensions, Hq, < © ^ 2ir with 

Dom(-He) = {ip G Dom(i^) : a_(V>) = -e ie a+(V>)}- (22) 

We take the self-adjoint operator Hq to be the realisation of the formal operator ( |10[ ). 

lr The notations / = O(g) and / <§C g both mean that there exists a positive constant C (the "implied 
constant") such that / ^ Cg. 
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3 Quasimodes 



3.1 Definitions and basic properties 

Let H be a self-adjoint operator in a Hilbert space, without continuous spectrum. 

Definition 3.1. A quasimode of H with discrepancy d is a pair (ip,^) £ Dom(ff) x R 
such that 

\\(H - ^\\ ^ dU\\. (23) 

We are interested in the situation when the quasi-eigenvalue and quasi-eigenfunction 
ip approximate true eigenvalues A,- and eigenfunctions <j)j of H. In this direction, the 
following classical results apply (see, e.g. [3"2"| 13]) 

For a quasimode with discrepancy d, the interval [fi — d, fj, + d] contains at least one 
eigenvalue of H. 

If we consider instead, the interval [n — M, /x + M] where M > 0, then 

E K^^>| 2 ^^ll^l| 2 . (24) 

In particular, if tp is normalised, and the interval [fi — M, fi + M\ contains only a single 
eigenvalue with eigenfunction <fi, then there is a phase x £ [0) 27r) such that 

(25) 

These results will be the main tools by which we relate the quasimodes constructed in 
the next subsection to the eigenfunctions and eigenvalues of H@. 



3.2 Quasimodes of delta perturbations 

Let / C 1 be a finite interval containing at least one point Ej of the spectrum of —A. Let 
a 6 [0, 1]. We will associate to the interval I a family of quasimodes parametrised by a. 

We first define 

^ ■= E + T^e Vlc ^ - e ' ' (26) 
where Vs is the spectral projection operator onto the set S, 

V S f:= E(/>^, (27) 

and I c is the complement to /. We have the following: 

Lemma 3.2. For z ^ Ej for any Ej € I, the function ip a ,i,z satisfies 

with the second term being bounded by a constant independent of I, z and a 6 [0, 1] . 
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Proof. We have 



E 



|3>i(p)i 5 



+ 



(7 



(^ ? --z) 2 |l-e i0 | 2 



(29) 



By (16) we get 



-i<9 



EJl-e 



-i©i 



#i(p) 
i(l + e" i0 : 



(1-e 



1 + 



*J-(P) 



sin© 1 



1 + £<? 1 - cos © 1 + #| 



#i(p). 



using 



i<9 



sin© 



1+e 

' 1 - e ie ~ 1 - cos ' 



By Parseval's identity 



sin© 



cos© J (1 + E]) 2 ' 



(30) 



(31) 



(32) 



Finally, to show that the right-hand side of (32) is finite and does not depend on /, we 
observe that it is bounded by 



f,/ sm© y \<P 3 ( P r 



i 



1+t 2 



sin © 
1 — cos 



dW(t), (33) 



writing the sum as a Riemann-Stieltjes integral. The spectral counting function N(t) was 

s, 



defined in (13). Integrating by parts, we get 

2 



sin© 



cos©; (1 + E]) 2 



If sin 
o dt \l + t 2 V l-cos© 



JV(t) dt. 



(34) 

Since iV(£) <C i by Weyl's law, we see that the integral in (34) is a finite constant. □ 
Let n = fi(a, I) be a solution to 



(7 



^(ffi-e'V.) (p). 



(35) 



Then the pair (ip a /, //) is a quasimode for where ^ j := ip a l a- This follows from the 
following proposition: 



Proposition 3.3. The function tp a j belongs to Dom(ffe) and satisfies 



||( J ffe-^ (J ,/|| 2 = (i- C T) 2 C/(0,/x) + ( 7 2 ^h+E j n+- 

Erfl V 



sin© 
— cos 



I^-(p)I 2 
(i + ^ 2 ) 2 ' 

(36) 
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Proof. First of all, let us prove that ip a i G Dom^©). 
We can write 



Ej- fj, 1 - e i61 



TMSi - e*Vi) + G» " e (37) 



Using ( 30 ) we can express this as 



= &,J + (<?i " e*Vi) , (38) 



where 



o" Ej <rsin(9 1 \ - . . , . . , , 

rri? + i-cos e i +g ; r» # 'W- (39) 



Now observe that due to the definition (35) of /i, tp a ,l(p) = 0, so ip^j G P p . Thus (38) 
justifies the assertion G Dom^©). 

Since -£Tq is an extension of -ff©, we have 

(H e - fj)^ = (Ho - + Y^I© (0 " + e i0 (i + . (40) 

Now, 



/ tt \ 7 sr^ I , aEAEj - a) a sin O Ej — u\ - - 



Y, (l + (W)£? + + ^(E, - ,,) (41) 



Using (16) we find 



<^,(i-/i) 5i + e ie (i + /i) 5 _i) 



again using (31). This leads to 



i + fi 



l + E] 



(42) 



1 — p~ ie / sin (9 

' •( 1 +^)-j4^(^-/i)j*i(p), 



^ (ife - - e^i) = -a £ ^1 + ^ + i^f^i " /*)) (43) 



3=1 



and combining this with ( |41[ ), we get 



E( _ sin 



gj(g) 

cos v J w i 1 + " J 



(44) 
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Since the summations in (44) are over disjoint sets it is easy to calculate the norm: 



ll(^e-M)^|| 2 = (l-<r) a X;i*i(P)l 5 



e&i V 



sin© 
— cos 



1 + E])f 



(45) 



□ 



3.2.1 Proof of theorem ll.il 



The first part of the theorem follows from lemma 3.2 and proposition 3.3 and the definition 
of a quasimode. 



For the second part, let fii ^ \X2 be two solutions of (35). We have 

T 2 



E 



sin© \ 2 \$j(p)\' 



cos©/ (l + E\ 



2\2' 



using (32). By elementary algebra, 

1 



(Ej - m)(Ej - fi 2 ) Hi - Hi \Ej - m Ej - H2 



so, since by (35 ) 



E 



E 



i*i(p)i s 



Ej - Hi £-? T Ej - H2 ' 



we get 



^(Ej-Hi^Ej-^) 



0. 



(46) 

(47) 

(48) 

(49) 
□ 



3.2.2 Controlling the discrepancy of quasimodes 

By tuning the parameter a, and choosing the interval I accordingly, we can find fix quasi- 
modes with particular properties. In the previous section we have seen that sets of quasi- 
modes with a = are orthogonal. We are particularly interested in when the discrepancy 



is small. In this subsubsection we prove corollary 1.2 that quasimodes with a = 1 can be 
made arbitrarily precise, and that quasimodes with a = also can have a simple bound 
for the discrepancy. 
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Proof of corollary Choosing I = [0, T] for T > E\ with a = 1 gives, by theorem 
that the discrepancy of ipij satisfies 



d 2 Uu\\ 2 



£j 1 + ^ + T^^-''> 



\Mp)\ 2 

(l+£ 2 ) 2 ' 



(50) 



By lemma 3.2 we see that the norm of ipij is bounded away from by a constant, so that 
the asymptotics for d are given by the term on the right-hand side of (50). Using Weyl's 
law, we can estimate 



sin 



\$j(p)\ 2 < v 



+ £ 2 ) 2 ^ T 



(51) 



which can be made arbitrarily small by increasing T. 
For the second part with a = 0, we have 



(52) 



We observe that splitting the sum in p5|) leads to 



£ I^'(P)! 2 



Ej£l 
Ej>ii 



Ej-n 



E 



i^-(p)i 2 

(i - Ej ' 



(53) 



Denote by EL and EL the largest and smallest points of the spectrum (Ej)JL 1 lying in 
the interval I. Then 



Ej>fi 



Ej>/i 



and 



£ |^-(p)| 2 ^ ( M -E_) E 



Ej</i 



Ej<» 



Adding these inequalities, and using (53), we get 

EiW<(^-^-)E§^ 

Ej£l E 3 £l j ^ 

Ej>fi 



fi-E, 



Ej<n 



(54) 



(55) 



(56) 



Since E + - E_ ^ we get 



2d 2 ||^o,/|| 2 < 



(57) 



Finally, 



^)E 



1^ " H 



w 2 E 



(58) 
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noting that \Ej — /i| $J 1(1) for \i 6 I. 

We now consider the case with a = 0, and / containing only the two levels Ej, Ej+\. 



We can solve (35) directly to get 



\<P J+1 (p)\ 2 E J + \<P 3 (p)\ 2 E J+1 
|^+i(p)| 2 + \<Pj(p)\ 2 

Substituting this value of [i into the definition of V'o,/ we get 
W - p + p -^i+i 

\<P 1 + 1 (p)\ 2 + \<P,(p)\ 2 f^_ _ + _J_ <[> 



^7+1 ~~ ^7 



i+1 



So, 



VII 2 = d2 IIVvl 



2I ^ +1 (p)| 2 + |^(p)| 2 



(Ej+i-Ej) 2 



1 1 

+ 



|^(p)| 2 l^+i(p)| 2 



4d 2 ||Vo,/|| 2 



using the fact that 
\<P H1 (p)\ 2 + \<P,(p)\ 2 



i(If 



_J_4_ 1 | -94- 



+ 



*i+i(p) 



(59) 



(60) 
(61) 

(62) 
(63) 

>4. (64) 

□ 



The existence of arbitrarily precise quasimodes can be used to give a new proof of 
the often-used representation for eigenvalues and eigenfunctions of rank-one perturbations 
(see e.g. [H EH E3 E] ) 



Theorem 3.4. The solutions A to the equation 



1 



E; 



e—x-YTe 2 )^ 2 



sin© ^ |^j(p)| 2 
1 - cos 1 4- E 2 ' 



are eigenvalues of H@ with corresponding eigenfunctions given by 



5=1 



(65) 



(66) 



Note that the left-hand side of (65) converges pointwise, and (66) converges in L (Q). 



By analysing the resolvent, it is possible to extend theorem 3.4 to get the following 
[281 Theoreme 2], 



Theorem 3.5. Apart from the solutions to (65), there are no other points of the spectrum 
of H@ in any of the intervals (Em, Em+i)- 
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4 Localisation results 



In this section we will consider the extent to which eigenfunctions of Hq can be approxi- 
mated by quasimodes. In particular we will focus on the quasimodes with a = 0. First we 



shall prove proposition 1.3, which is straightforward. Then we shall show that strength- 
ening the assumptions made on the spectrum of A leads to a proof of full convergence. 

Proof of proposition [7T3| The length of the interval / is £(I) = E c —E^. Let M = minj^ - 
jjb, fj, — Eg} ^ minjE'd — E c , E^ — E a }. By applying (24) with this M we get, 



y 0,/| 



A 3 -£[B a ,B d ] 



Amm{E d -E c ,E h -E a yj ' 



(67) 



From theorem 3.5 there are only three eigenvalues of H@ in the interval [E a , Ed]. It there- 
fore follows that for at least one of these three eigenfunctions its inner-product squared 
with i is at least i of the right-hand side of ( 67 ) . □ 



We now consider how to improve proposition |1.3| at the expense of making further 
assumptions about the spectrum of —A. For simplicity we will focus henceforth on the 
choice of parameter = ir. 



Eg. 



E h 



A 



E r 



Figure 1: Part of the spectrum of —A and H n . Vertical bars denote eigenvalues of H n , 
and circles denote eigenvalues of —A. The triangle fj, is a quasi-eigenvalue approximating 
A. See main text for further explanation. 



In figure[T]a cartoon of part of the spectrum of H n and — A is displayed. Highlighted are 
four consecutive eigenvalues of — A, labelled E a ,E^, E c and Ed, chosen so that E c — E^ ^ e. 
(The positions of all points depend one.) 

Between E\> and E c is an eigenvalue, A, of H n . 

We find a quasimode ipoi associated to the interval I = [E^, E c ] with quasi-eigenvalue 
\i approximating A. By corollary 1.2 the discrepancy of this quasimode is no greater than 
e/2. Between E c and Ed is another eigenvalue A* of H n . In order to be able to apply (25), 
we need to be sure that A* is not too close to E c . An argument to show that this is the 
case is given below. 

The eigenvalue between E a and E^ can be handled with a similar method. 

We shall make the following assumption on the spectral sequence of A: 
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Assumption 4.1. For some < q < 1/2 and 1 < p < 2(1 — q), there exists a sequence 
{e n ) ( ^L l , £„ | such that for each n there are four consecutive eigenvalues, E^n) < 
£"b(n) < E c (n) < E&{n) <C £n P , satisfying 



E c — E\> -C e n 
£ d - E c » 4 



(68) 



as n 



oo. 



Assumption 4.1 asserts that the positions of eigenvalues of —A occur with the spacings 
as described above, and furthermore that this does not happen too high up in the spectrum. 
This upper bound is necessary as a consequence of the non-uniform convergence in A of the 
series in ( |65[ ). In appendix B we show that assumption 4.1 is satisfied almost surely if the 
sequence (Ej) comes from a Poisson process. In this sense, assumption 4.1 is consistent 
with the Berry- Tabor conjecture, if —A is the Hamiltonian corresponding to an integrable 
dynamical system. 

We shall also assume a lower bound for the absolute values of the eigenfunctions <&j at 
the point p. 

Assumption 4.2. There exists a constant cq > independent of j such that 



l*i(p)| > co- 



(69) 



Remark 4.3. In fact we require only that assumption 4-2 holds for (possibly a subsequence 



of) the sequence of pairs 'Pbip) and &c(p) f° r eigenfunctions associated with the sequences 



of energy levels E^ and E c defined in assumption 4-1 



We recall that the spectral sequence is defined in such a way that $j(j>) ^ for all 
j. Thus assumption 4.2 disqualifies subsequences of eigenfunctions converging to at the 
point p. 

This assumption reflects the fact that if becomes small, two eigenvalues of H n 

will approach Ej. Then we would only be able to prove that the quasi-eigenfunction 
approximates a certain linear combination of these eigenfunctions of —A, rather than an 
actual eigenfunction. Assumption 4.2 can be relaxed slightly — see remark |4.5| below. 



Theorem 4.4. Assume that the spectrum of —A satisfies assumptions 4-1 and U^M Then 
the sequence of quasimodes tftoj associated to the sequence of intervals I = [E^,E C ] and 



[i € /, with E a , . . . , E& as described in assumption \4-l\ after normalisation, converge in 
I? to a subsequence of true eigenfunctions of H n . 



Let us fix a point n of the sequence (e n ) with e n 
statement of theorem 14.41 



e, and I fixed as described in the 



Proof of theorem 4-4- I n order to use (25) we will employ partial summation, to estimate 



the position of eigenvalues of H n . If g is a smooth function, then 



£ g{E 3 )\^{p)\ 2 = g(Y)N(Y) - g(X)N(X) 



Y 



X 



g'(t)N(t)dt, 



(70) 
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where N(t) has been defined in (13). Equation (70) may be proved by Riemann-Stieltjes 

I b 

t 



integration. Let A* be the solution of (65) lying between E c and E d . Let 

1 + tA* 



g(t) ■■- 



1 



t-X* 1 + t 2 (t- A*)(l + t 2 )' 
and observe that g(t) > if t > A* and g(t) < if t < A*. By d65j) we have 



= y £g(E j )\$ J (p)\ 2 ^g(E c )\<P c (p)\ 2 + £ g(E 3 )\<I> 3 (p)\\ 

3=1 



Now, by ([TOp , 

^ 5 (^)|^(p)| 2 = -g(E d )N(E d ) 



Ej^Ed 



g'(t)N(t)dt 



E d 



— jf^ g(t)dt + 0i K g(E d )E 1 d /2 + J e \g'(t)\t l l 2 dt 



using (13). 
Since 



we get 



-1 



1 



+ 



2t 2 



\g'(t)\tV 2 dt 



(t-X*) 2 1 + t 2 (1 + t 2 ) 2 ' 

t l/2 



■a (*"*) : 



dt 



£ d -A* 



1 



(t - A*) 3 /2 



dt 



E d 



1/2 



i^d-AV (£ d -A*)V 2 



1/2 



£ d -A* 



We can also calculate 

poo 

/ s(t)dt 

So we have 



E, t-X* 1 + t 2 



dt 



In 



x; ^)i*i(p)i s 



1, E d -X* 



47T 



In 



+ 



E d -X* 



(71) 



(72) 



(73) 



(74) 



(75) 



(76) 



(77) 



in which the dominant term on the RHS is actually the error term. We have, from ( 72 ) 

,1/2 



F 

g(E c )\<P c (p)\ 2 ^ 



(78) 



£ 



-p/2 



I*C(P)| ! 

E c - X* ~ E d -E c - (X* - E c 



ei - (X* - E c ) ' 



(79) 
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implying the lower bound 

A* - E c > e p/2+q . 

To see this, observe that if A* = o(e q+p ^ 2 ) then we would have from (79) 



\Mp)\ 2 

E c - A* 



-p/2-q 



(80) 



(81) 



a contradiction. 



By the same method, we can e stablish the same bound for the solution to (|65|) between 

p/2 



3.5 



we deduce that there is an interval of size M x e q 



E a and E^, and by theorem 
about fi such that [/j — M, \i + M] contains only one eigenvalue of H w . Since q + p/2 < 1, 



and since the discrepancy of ipoj is 0(e), equation (25) allows us to conclude that the 
normalised quasimode differs from the true eigenfunction associated to \x by an amount 
which converges to as e — > 0. □ 



Remark 4.5. From the proof of theorem 4-4 we see that we can relax assumption 4-2 to 
demanding only that \<Pj(p)\ 3> e r l 2 with 0<r<l — q — p/2. However, in a generic situ- 
ation this is unlikely to be achieved. In appendix\A\ we show that for a badly- approximable 
position of the point p in a rectangle, the best possible bound is 



\$j{p)\ » 



1 



(82) 



which is not sufficiently slow. 



5 Application to rectangular Seba billiards 



In this section we will apply theorem 4.4 to the original Seba billiard [12]. We consider 
a rectangular billiard Q = (0, 2a) x (0, 26) C M 2 and point p = (a, b) at the centre of the 
billiard. However, we remark that we could position p at any point with co-ordinates that 
are rational multiples of the side lengths without significant changes to the forthcoming 
analysis. 

The eigenvalues of —A, the Laplacian with Dirichlet boundary conditions, are given 

by 

7r 2 / n 2 m 2 



^• m -Tl^ + ^ 2 " ) ' (83) 



where n, m G N, and the corresponding eigenfunctions are 

. . 1 /nirx\ / 771717/ \ , .s 

^ m (x, y ) = - 7 =sm( — )sm( — ). (84) 

If either n or m are even, then the symmetry of the problem forces $ n m{p) = 0. So for 
these values of n and m, <& n ,m £ D p , and are automatically eigenfunctions of the extended 
operator H n . We exclude these eigenvalues from the spectrum, as discussed in section [2} 

Instead, we concentrate on the more interesting subsequence where n and m are both 
odd, e.g. n = 2s + 1 and m = 2t + 1 with s, t = 0, 1, 2, . . . Then we have 

<2> g) ,(p) = _L(-1)*+*, (85) 
\J ab 
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so that along this sequence assumption 4.2 is satisfied. The corresponding set of eigenvalues 
is given by 



E s .t = ^ 



+ 



b 2 



s,t = 0, 1,2, 



(86) 



For generic choices of a and b it is conjectured that the set of values given by (86) behave 
statistically like the event times of a Poisson process |17} 134"! [351 136] . Under assumption 



4.1 for the set of values (86), theorem 4.4 asserts the existence of a subsequence (j n ) C N 

Il0in _ V'nll ->• as n oo, (87) 



such that 



where </>,■ are eigenfunctions of H^, and ifi n are of the form 



i>n = j= (% + (-l)^^V 



+1 > 



where /3 n can be or 1 and depends on the relative signs of $j n (p) and $j n+ i(p). So the 
subsequence (4>j n ) converges to a superposition of two consecutive unperturbed eigenfunc- 
tions of —A. 

The consequences for this subsequence are most striking when one considers the mo- 
mentum representation. This is is given by the Fourier transform; 



1 

2^ 



(89) 



— OOJ — OO 



For an ergodic system, the quantum ergodicity theorem of Snirel'man, Zelditch and Colin 
de Verdiere |37} 155} [39] would imply that the momentum representation of almost all 
eigenfunctions equi-distributes around the circle of radius y/Ej as j — > oo; 



1 



(x)\ — > — S(\x\ — Ej) as j — > oo, 



TT 



(90) 



where convergence in (90) is in the weak sense. Seba billiards are not ergodic, but we see 
a very different behaviour to (90) for the subsequence (4>j n )- 

From Parseval's theorem, it follows that 

4>j n — ip n — > in 1? norm. (91) 
The momentum representation of the unperturbed eigenfunctions <S> n ,m is 



&n,m(Px,Py) 



Since 



2irnmv ab 



(4:p 2 a 2 — n 2 ir 2 )(4:p 2 b 2 — m 2 Tr 2 ) 



-jr ((-l) n e- 2iap * - l) ((-i) m e - 2ibp y - 1 



nir 



1 



1 



1 



<ip x a 2 — n 2 Tr 2 2 \2p x a — nir 2p x a + nn 



(92) 
(93) 



we re-scale and write 



@n,m(np x ,mp y ) 



7rva6 
2nm 



(5 n (2p x a + 7r) - 5 n (2p x a - tt)) (S m (2p y b - it) - 5 m (2p y b + tt)) 



(94) 
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where 5 n is the smoothed-delta function 

1 _ p-int 

Sn(t)--= — • (95) 

7T1C 

The function 5 n (t) converges weakly to 5(t) as n — > oo. Furthermore, it satisfies 

2n 

IM*)| 2r 8(t) asn-^oo. (96) 

TT 

Hence 

nm\$ n , m (np x , mp y )\ 2 ~ ab (5(2p x a - it) + <5(2p x a + 7r)) (5(2p y b - tt) + <5(2p y 6 + 7r)) (97) 
as n,m -> oo. The momentum eigenfunction localises around the 4 points 

which satisfy p%+Py = E n ^ m . Since ip n is a superposition of #j n and ^j n +i, the states 
in the subsequence 0j n become localised around 8 points, which all lie on the circle with 
radius \/~Ej, very much in contrast to the expected equi-distribution (90) for ergodic 



systems. Numerical simulations illustrating this behaviour have been presented in [26 . 
This localisation is, in some sense, analogous to the scarring phenomenon which occurs 
in some chaotic systems. Since these states are not associated with an unstable periodic 
orbit, they do not fall into the very precise definition of a scar given in |40j . Rather 
they are localising around ghosts of departed tori of the unperturbed integrable system. 
Nevertheless they cannot be explained simply by using torus quantisation, and so they 
provide a further example of the already rich behaviours in systems with intermediate 
statistics. 
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A Non-constant unperturbed eigenfunctions at the position 
of the scatterer 

In order to consider what can happen when the value of the unperturbed eigenfunctions 
at the position of the scatterer can vary, let us consider the rectangular billiard 17, with 
sides of length a and b, and Dirichlet boundary conditions. 
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The energy levels are given by 

E = E n , m = TT 2 ( ^2+^ ) • (99) 



2 2 

n m 



for n, m 1 integers. 
Lemma A.l. 



1 4vr 4 1 . , 

> 7212^2- (100) 



n 2 m 2 a 2 b 2 E 2 



Proof. We have 



o 2\2 44 22 44 

, n m \ n it n m A m tv 



a 2 b 2 J a 4 a 2 b 2 6 4 



4 4 
n 7T 


2 2 

n m 


a 4 


2 a 2 b 2 




2 2 


= £ 2 - 


. n m a 

4 7T 


a 2 6 2 



n— — r 



(101) 

and then re-arrange to get the required estimate. □ 
The eigenfunctions themselves are proportional to 

* (=?)-(=?)• < 102 > 

Let us choose the point p = (x p ,y p ) £ Q so that x p /a and y p /b are badly-approximable, 
in the sense that 

C 

Vn,rGZ, (103) 

n 

(this is the best we can hope to do if we want to bound the eigenfunctions away from 0). 
Then 

n^=r + ??(ra) (104) 
a 

where •& can depend on x p and r and satisfies 

|0(n)| > - (105) 
n 

uniformly. Furthermore this bound is achieved if r/n is a continued fraction approximant 
to Xp/a. We get 

sin 2 ( -)»^. 106 

With a similar bound for the contribution of the y-coordinate, we find that the best bound 
we can obtain is 

\&n,m{p)\ 2 » > ~B2 (1° 7 ) 

n A m z E z 

and this bound is sharp. 



B Assumption 4.1 for the event times of a Poisson process 



The purpose of this appendix is to prove the following result. Let < q < 1/2 and 
1 < p < 2(1 — q) be fixed throughout. 
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Proposition B.l. Let P = {Ej)'^ =l be the sequence of event times for a Poisson process 
with parameter 1. There is, almost surely, a sequence (e n )^ =1 , e n I such that for each n 
there are four consecutive members of P , < E^ < E c < E d < e n p , satisfying 

E c — Efr < e n 

E d -E c > el (108) 
E h -E a > el 



Thus, assumption |4.1| is almost surely satisfied for a Poisson process. 

As a model for a Poisson process, we will let 6>6> ... be a sequence of independent 
exponentially distributed random variables with parameter 1. Then, defining 

£1 = 6 

£2 = 6+6 (109) 



The sequence P = (Ej) c *L 1 so- formed is a Poisson process. 

Proposition B.2. Let e > 0. The probability that there are four consecutive members of 
P, E & < E h < E c < E d < e~ p satisfying 

E c -E h <e 

E d -E c > e q (110) 
E h -E & > e q , 

is 1 - 0(e°°). 

The notation 0(6°°) refers to a quantity which goes to zero faster than any power of e. One 



can say that the event described in proposition B.2 occurs with overwhelming probability. 

Let us fix 1 < p' < p, and chose N = 3M ~ e~ p where M £ N. Let us define the 
events Sj, j = 0, . . . , M — 1, by 

Sj = {6i+i > e\ 6i+2 < e, 6i+3 > e q }- (HI) 



Lemma B.3. The events Sj , j = 0, . . . , M — 1 are independent, and the probability that 
at least one of them occurs is 1 — 0(e c 



Proof. The independence of the events Sj clearly follows because they are defined on 
independent random variables. We first calculate the probability of one of them. By 
independence of 6>6>6> 



F(S ) = P(6 < £)F(6 > e 9 )IP(6 > e 



\ / roc \ 2 

e~ x dx I [ / e~ x dx 



/ \Jei 



(l-e-)(e- 9 ) 2 

e + 0(e l+q ). (112) 
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Then, by independence of the SjS, 

pi := P(at least one Sj occurs) = 1 - (1 - P(5 )) M 

= 1 - (1 - e + 0(e 1+q )) M . (113) 

So, we have 

log(l - pi) = -Me + 0(Me 1+q ) ~ - Je 1 ^'- (114) 
For e sufficiently small, this yields 

1 - pi < exp (-1 £ -(p'-i)) = 0(e°°). (115) 

□ 

The probability that the upper bound of e~ p is met is given in the following lemma 
Lemma B.4. The probability that En < e~ p is 1 - 0(e°°). 

Proof. Let a > 0. The probability density for E N is r(N)- 1 x N - 1 e~ x . So 

1 f'°° 

P2 := P(^iv < N 1+a ) = 1 - — — / x N ~ l e~ x dx 

= i - exp( ~^! +Q) rix+N^N-^-'dx 

Jo 



r(N) 



N-l 



1 ~ exp(-iV'-)Arf' + "»(^) £ r(w _ j ' )Ar , 1+ „, J . (116) 



1 AT-? 

expanding the binomial. Using ^ we can estimate 

r(N-j) r(N) 



N-l , JV-1 



V 1 1 V ^ 1 1 

r(N - j)N( 1+a )i ^ r(N) < r(N) ' ^ ^ 

where the implied constant could depend on a. This leads to 

exp(-iV 1 + Q )iV( 1 +«)(^ 1 ) 

7^) 

exp(-iV 1+Q + N)N a ( N ~ 1 ') 
y/2ir(N - 1) 

« (118) 

where Stirling's formula has been used. This last line is 0(e°°) since N- 1 ~ e' 5 '. Finally 
setting 

a = 4~l>0 (119) 
P 

gives the required estimate. □ 



20 



Proof of proposition \E~2 , We are interested in the events corresponding to lemmata B.3 



and B.4| happening simultaneously. By the inclusion-exclusion principle, the probability 



that this happens is at least p\ + p2 — 1 = 1 — 0(e°°). □ 



Proof of proposition B.l. Let S n , n G N be the event that there are found four consecutive 
members of P, E a < E^ < E c < < n p satisfying 

E c — E^ < — 
n 

E d -E c >^- q (120) 
^b-£a>-- 



By proposition B.2, P(£^) -C ~i- Hence, by the Borel-Cantelli lemma, the probability that 
infinitely many £^ occur is zero. Equivalently, only finitely many £^ occur, almost surely. 
So, almost surely, there is an infinite subsequence of n G N such that £ n occurs. □ 
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